function[c,f,s]=pdefun(x,t,u,ux)

c = [1;1];

p = 1+0.5*abs(sin(pi*t));

y(1)=u(1)*(1.8+0.3*abs(sin(pi*t))-(0.3*abs(sin(pi*t)))*u(2)*inv(1+0.4*abs(sin(pi*t))*u(1))-u(1));

y(2)=u(2)*(1+0.5*abs(sin(pi*t)))+(0.2*abs(sin(pi*t))*u(1)*inv(1+0.4*abs(sin(pi*t))*u(1))-u(2));

F=[y(1);y(2)];

s=F.*[1;1];

f = [0.14*inv(p^2)*ux(1);0.10*inv(p^2)*ux(2)];

end